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ABSTRACT 

Observations have revealed prodigious amounts of star formation in starburst 
galaxies as traced by dust and molecular emission, even at large redshifts. Recent 
work shows that for both nearby spiral galaxies and distant starbursts, the global 
star formation rate, as indicated by the infrared luminosity, has a tight and almost 
linear correlation with the amount of dense gas as traced by the luminosity of 
HCN. Our surveys of Galactic dense cores in HON 1—0 emission show that this 
correlation continues to a much smaller scale, with nearly the same ratio of 
infrared luminosity to HCN luminosity found over 7-8 orders of magnitude in 
Ljn, with a lower cutoff around lO'*'^ Lq of infrared luminosity. The linear 
correlation suggests that we may understand distant star formation in terms of 
the known properties of local star-forming regions. Both the correlation and the 
luminosity cutoff can be explained if the basic unit of star formation in galaxies 
is a dense core, similar to those studied in our Galaxy. 
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1. Introduction 

Recent work has revealed large amounts of dust and molecules in starburst galaxies, 
even at large redshifts (e.g. Isaak et al. 2002; Reuland et al. 2003; Greve et al. 2005; 
Solomon & Vanden Bout 2005). Understanding star formation in galaxies at high redshift 
is a critical step in understanding the formation of galaxies in the early Universe. 

The simplest and most widely used relations between the star formation rate and a 
property of the interstellar medium (ISM) are the so-called "Schmidt laws". Schmidt (1959) 
proposed that the star formation rate was proportional to p^, where p is the gas volume 
density. In their modern form, these "laws" relate the surface density of star formation to 
the surface density of gas: 

^SFR = AEj, (1) 

(e.g., Kennicutt 1998). The index N has been inferred by various authors to be in the range of 
1 to 2 (Kennicutt 1997). Measurements of HI, CO, and Ha on a large sample of normal spiral 
galaxies and starburst galaxies (Kennicutt 1998) could be fitted with N — 1.4. It should be 
noted that for many of these galaxies, particularly the starburst galaxies, in the Kennicutt 
sample the surface density or surface brightness was not measured but obtained only by 
dividing the luminosity by a characteristic size often obtained from another parameter. 

On a global scale, including luminous and ultraluminous starburst galaxies, there is 
a correlation between the total luminosity of far-infrared emission, which traces the star 
formation rate (e.g. Sanders and Mirabel 1996; Kewley et al. 2002), and the total luminosity 
of CO, tracing the molecular gas mass. However, this relation is not linear; the ratio of L/r 
to Leo increases with increasing L/^ ( Sanders and Mirabel 1996; Solomon et al. 1997; Gao 
& Solomon 2004a, 2004b; Solomon & Vanden Bout 2005). Does this variation reflect an 
increase in "efficiency" of star formation, or an increasing failure of CO to trace the gas that 
is relevant to star formation? 

The latter possibility is suggested by data on HCN J = 1 — in galaxies. A recent 
survey of HCN J = 1 — emission in 65 normal spiral and starburst galaxies found that the 
star formation rate, as traced by the infrared luminosity, has a tight and linear correlation 
with the luminosity of HCN (Gao & Solomon 2004a, 2004b). Those authors argued that CO 
is not hnearly correlated with star formation because it traces only the low density CMC 
envelopes, not the really active star-forming part, the dense cores. The critical molecular 
parameter that measures star formation rates in galaxies is the amount of dense molecular 
gas, measured by the HCN luminosity. Because HCN traces the dense gas better than CO, 
^HCN has a tighter correlation with Lm than does Leo- The correlation remains linear over 
a factor of 10^ in luminosity for both normal galaxies and extreme starbursts, like luminous 
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and ultraluminous infrared galaxies (LIRGs and ULIRGs, see Sanders and Mirabel, 1996). 
Gao and Solomon (2004a) therefore argue that both normal galaxies and starbursts should 
have the same star formation rate per amount of dense gas. 

At lower luminosities, the relation for galaxies between Lir and Leo is linear, but the 
ratio between L/r and Lqo has a dispersion of an order of magnitude. An even larger 
variation in Lir/Lqo (several orders of magnitude) is seen in Galactic clouds (Mooney & 
Solomon 1988; Evans 1991; Mead et al. 1990). The dispersion in the relation for galaxies is 
less if Lhcn is used instead of CO; the same is true for other tracers of dense gas in Galactic 
cores (e.g. Mueller ct al. 2002, Shirley ct al. 2003), suggesting that studies of Galactic star 
formation can shed light on the trends in other galaxies. 

In the Milky Way, star formation is dominated by clustered star formation (Carpenter 
2000; Lada & Lada 2003 ). Clustered star formation produces stars with a range of masses, 
but massive stars form nearly exclusively in clusters within massive dense cores. Since 
massive stars dominate the luminosity, they are the stars directly probed in studying star 
formation in other galaxies. Thus, understanding the relation between star formation and 
massive, dense cores in our Galaxy may shed light on the starburst phenomenon. 

These cores are dense, turbulent, and dusty. They are well identified by the continuum 
emission from dust (e.g., Mueller et al. 2002) and line emission from molecular dense gas 
tracers hke CS (Plume et al. 1992, 1997; Shirley et al. 2003) and HON (Wu & Evans 2003). 
In our previous work on a large sample of massive cores, we found that the bolometric (almost 
all far-infrared) luminosity of the cores is roughly proportional to the mass inferred from the 
dust emission (Mueller et al. 2002) and the virial mass determined from CS (Shirley et al. 
2003). This result suggests that a relation between Lm and Lhcn may exist in Galactic 
dense cores, possibly similar to that in starburst galaxies. 

One difficulty is that systematic studies of the Galactic dense cores have used dust 
continuum emission and CS, while HCN has been more commonly used for studies of galaxies. 
To facilitate comparison with the HCN J = 1 — galaxy survey, we have mapped the HCN 
J = 1 — transition in a sample of 47 Galactic star forming cores (Wu et al. 2005, in 
preparation). In this paper we will summarize the results from this survey in connection 
with the HCN surveys of other galaxies. 

The sample mapped in HCN J = 1 — is mostly a subset of a larger sample of Galactic 
massive dense cores that have been mapped in CS and dust emission (Shirley ct al. 2003; 
Mueller et al. 2002). The sources in this category have infrared luminosities ranging from 
10^ - 10^ L© and most contain compact or ultracompact H H (UCHH) regions. To extend 
the sample towards lower luminosities, we selected 14 IRAS sources from outflow surveys 
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(Zhang et al. 2005; Wu et al. 2004) and a few from other publications. 

2. Observation &; Data Analysis 

Observations of HCN J = 1 — > (88.6318473 GHz) on Galactic dense cores were made 
with the 14-m telescope of the Five College Radio Astronomy Observatory (FCRAO) in 2004 
April, December and 2005 February. The FWHM of the beam for this frequency is 58". The 
16-element focal plane array (SEQUOIA) was used, with typical system temperatures 100- 
200K. A velocity resolution of 0.1 km s"^ was achieved with the 25 MHz bandwidth on the 
dual channel correlator (DCC). We convert the measured to Tr via Tr = T\l {rjpssVc)-i 
with r]pss = 0.7. The value of ric depends on source size; for the typical map in this study 
(~ 10'), ric = 0.7. The map size was extended until the edge of the HCN J = 1 ^ emission 
was reached, typically at the 2a level (mean a ~ 0.3 K km/s), so we could get the total 
HCN luminosity. 

Maps of HCN 1-0 emission were obtained for 47 Galactic star-forming cores. Detailed 
results will be presented elsewhere (Wu et al. 2005, in preparation), but we give a few 
properties of the sample here, which are relevant to this paper. More than 90% of the cores 
were well resolved by the maps. We focus on these resolved sources in this paper. 

The size of the core is characterized by the nominal core radius after beam deconvolution, 
Rhcn-i the radius of a circle that has the same area as the half peak intensity contour: 

A , fl2 

rlHCN = Jy[ — > , 

TT 4 

where A1/2 is the area within the contour of half peak intensity. The median Rhcni-o of 
the sample is 0.48 pc. 

The HCN line luminosity of each core, assuming a gaussian brightness distribution for 
the source and a gaussian beam, is: 

Lhcn = 23.5 X 10-6 xD'x (^^) x (^?-t^) x j T^dv (2) 

Here D is the distance in the unit of kpc, 6s and Obeam are the angular size of the source and 
beam in arcsecond. This method is parallel to that of Gao & Solomon (2004b), but adapted 
to Galactic cores. Lhcn ranges from 0.4 to 8000 K km s^^ pc^, with the median value of 80 
K km s""^ pc^ . 

The total infrared luminosity (8-1000 /xm) was calculated based on the 4 IRAS bands 
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(Sanders and Mirabel 1996), as was done for the galaxy sample of Gao & Solomon (2004a): 

LiR = 0.56 X X (13.48 x /12 + 5.16 x /ss + 2.58 x /go + /loo), (3) 

where is the flux in band x from the four IRAS bands in the units of Jy, D in kpc, and 
LiR in Lq. 

3. Compcirison of Milky Way and galactic Relations 

The derived Ljr and Lhcn are plotted in a log-log diagram in Fig. 1 to compare 
with data on galaxies from Gao & Solomon (2004a). The correlation of Ljr-Lhcn extends 
from galaxy scales to the much smaller scales of Galactic molecular cores. The fit for Gao 
& Solomon's galaxy sample is \og{LjR) = 1.00 x log^LncN) + 2.9, without a few galaxies 
that only have upper limits to HCN emission. This linear correlation continues to the 
Galactic massive cores, but a decline in Ljr occurs at 10^'^ Lq, below which the slope of 
the correlation becomes much steeper. When fit to Galactic cores with Ljr > 10"^'^ Lq, 
the linear least squares fit gives log{LiR) = 1.02(±0.06) x log(LHC7v) + 2.79(±0.16). This 
relation agrees remarkably well with the relation for galaxies, as seen in fig. la. In fig. lb, 
we fit simultaneously the data of galaxies from Gao & Solomon (2004a) and Galactic cores 
with Ljr > 10^'^ Lq. The result is log{LiR) — 1.01 x \og{LHCN) + 2.83, with a correlation 
coefficient of 0.99. Understanding the physics behind this linear correlation will lead to a 
better understanding of star formation on galactic scales. 

The linear \og{LiR) — \og{LHCN) correlation and the turnoff can be seen more clearly 
from fig. 2, where the distance independent ratio Ljr/ Lhcn has been plotted versus Ljr 
(top) and against Lmoi (bottom). A constant mean value of Ljr/ Lhcn is seen over 8 or- 
ders of magnitude in Ljr, from galaxies to Galactic cores, as long as Ljr > L{min), with 
L{min) being aroimd 10"^'^ Lq. The corresponding cut-off value for Lhcn is Lunit = 10^'^ K 
km s^"*^ pc^. For Galactic cores, {Ljr/ Lhcn) = 911 ± 227, with median 380. For galaxies, 
{Ljr/ Lhcn) = 950 ± 76, with median 855. The uncertainties are the standard deviation of 
the means, which are remarkably similar, though the dispersion is much higher for Galactic 
cores. The median for dense cores is significantly smaller than its mean, which indicates that 
the mean is dominated by a few quite large values. The logarithmic mean and median are 
2.64±0.53 and 2.54±0.56 for dense cores with Lhcn > Kniu and 2.91±0.24 and 2.93±0.25 
for galaxies. 

For comparison, we added CO data on Galactic cores (Mooney & Solomon 1988), galax- 
ies (Gao & Solomon 2004a) and high-z molecular emission fine galaxies (Solomon & Vanden 
Bout 2005) in fig. 2. The ratio, Ljr/ Leo, changes by two orders of magnitude as Leo 
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increases from Galactic cores to distant galaxies, confirming Gao & Solomon's conclusion 
that CO is not as good a tracer of star- forming gas as is HCN, especially for very luminous 
star burst galaxies. 

4. Discussion 

The fact that Ljr/Lhcn is similar, on average, for star-forming cores in the Galaxy, 
normal spirals, starbursts, and ULIRGs suggests the possibility of interpreting intense high-z 
star formation in terms of nearby high mass star forming regions. Before we can exploit this 
possibility, we must understand some key points. Why does Ljr/Lhcn rise steeply with 
Lrcn and then remain constant for Lhcn > Lunit'- And why is the ratio, Lir/Lco, NOT 
constant for starbursts? 

As a first step, we seek a more physical basis for the relations. We have so far discussed 
Lhcn as a measure of the mass of dense gas, but can we quantify this assumption? A roughly 
linear correlation between the mass of dense gas and bolometric luminosity has been found by 
our work on CS and dust emission (Shirley et al. 2003; Mueller et al. 2002). To see whether 
this applies also to HCN J = 1 — 0, we have calculated the virial mass of the dense gas 
{M {dense) ) using the most optically thin line (C^^S J = 5 ^ 4) to measure the linewidth and 
compared M{dense) to Lhcn- Based on the 31 cores with available C^^S J = 5 — > 4 data, we 
obtained the correlation: log(M^j^(i?i/cJvi-o))=(0.81±0.03) x log(Li/cJvi-o) +1.29(±0.09). 
The correlation is roughly linear. The Mmr-LncNi-o plot and details of the observations 
of C^'^S J = 5 — > 4 will be presented by Wu et al. (2005, in preparation). Leaving out one 
peculiar source (G35.58-0.03), we get {M {dense) / Lhcn) = 11 ±2 Mq/K km s~"^ pc^, where 
the uncertainty is the standard deviation of the mean; the median value is 6 Mq/K km s~^ 
pc^, indicating that the mean is affected by some quite large values. The logarithmic mean 
is 7±2 Mq/K km s~^ pc^. Some of the scatter in the ratio may be caused by distance 
uncertainties because the virial mass depends linearly on distance, while Lhcn oc D^. 

Even after establishing that Lhcn traces the mass of dense gas, it is not at all clear 
why Lin/ Lhcn should be constant, since the luminosity of a cluster is typically dominated 
by its most massive members, and L^, oc M^", with a ^ 3 to A. Indeed, below the cutoff of 
10^-^ Lq, Ljji does rise rapidly with M {dense). What causes the transition to a constant 
value? 

To solve this puzzle, we propose the existence of a basic unit of cluster formation. For 
M{dense) less than the mass of this unit, Lm/M {dense) rises rapidly with M {dense), as 
higher mass stars can form. For M{dense) greater than the mass of this unit, the IMF is 
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reasonably sampled and further increases in mass produce more units, but no further change 
in Ljr/ M {dense). If we suppose that larger scale cluster formation is built up by adding 
more and more such units, then the linear correlation between the total L/^ and M{dense) 
is a natural result. In that case, the only difference between star formation on different 
scales and in different environments-big clusters, normal galaxies, massive ULIRGs-is just 
how many such cores they contain. This is of course somewhat of a simplification because 
the cores have a range of properties. The fact that the mean is roughly twice the median 
for both Ljji/ Lhcn and M{dense) / Lhcn for Galactic cores suggests that the most extreme 
cases may dominate when averaged over a whole galaxy. Since most of our cores contain 
compact H n or UGH 11 regions, which trace the most massive stars, it is interesting to 
compare to the luminosity function of these regions. A study of Galactic UGH 11 regions did 
find a peak luminosity (~1O^L0) that is close to that of a basic unit (Cassassus et al. 2000). 
The detailed discussion of this model will be presented in a separate paper (Wu et al. 2005, 
in preparation). 

Another question is why docs Lj^jLco increase as we move from Galactic Gores or 
normal spiral galaxies to starburst galaxies. In Galactic clouds, GO can be used as a tracer 
of the overall mass of molecular clouds, even though it is optically thick and thermalized. 
However, GO does not specifically trace the mass of dense cores. For that, dust continuum 
emission, or molecules, like HGN, that are only excited at higher densities are required. In 
the Galaxy, these dense cores occupy a small fraction of the mass of a cloud (typically a 
few percent). If the fraction of dense gas in the overall cloud stays roughly constant, the 
relation between Ljji and Lqo can stay linear, but this relation is secondary to that between 
L/R and M (dense). For starburst galaxies, the fraction of the molecular gas concentrated in 
dense cores increases, causing (Gao & Solomon 2004a) the secondary relation between Ljr 
and Leo to become non-linear. 

These considerations lead us to offer some new versions of "Schmidt Laws" : 

M^MoT/r"^) - 1.4 X 10-^LHCN{Kkm s'^ pc^) (4) 

M^(MQyr-^) ~ 1.2 x 10'^ M (dense) (Mq) (5) 

where M^, is the star formation rate, and we have assumed that M^(MQyr~^) — 2.0 x 10~^^ Liji(Lq) 
(see Gao & Solomon (2004a)) and we use the fit to both dense cores and galaxies. The co- 
efficients (1.4 X 10~^ and 1.2 x 10~^) are very similar to, but shghtly less than, those values 
(1.8 X 10"''' and 1.8 x 10~^) given by Gao & Solomon (2004a), based only on the galaxies. 

Given these simpler, linear relations, how do we understand the usual star formation law 
(Kennicutt 1998) relations, with a steeper dependence (T,sfr oc T,^gas) of M* on gas mass? 
We suggest that the steeper dependence reflects the dependence of M(dense) on the surface 
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density, or total mass, of gas. All the tracers of gas used by Kennicutt trace lower density 
components, not the actual gas that is directly involved in star formation. Once it is clear 
that it is the dense gas mass that indicates the star formation rate, it becomes clear why 
the total surface density of gas may not be a clean star formation indicator. For example, 
there is no evidence that HI emission in galaxies correlates at all with the star formation 
rate deduced form the far-infrared luminosity, so its contribution to the surface density may 
have no effect on the star formation rate. Attempts to provide a theoretical framework (e.g., 
Li et al. 2005, Krumholz & McKee 2005) for the Kennicutt relations should also be able to 
explain the new relations. 

It will be interesting to sec how far these ideas can be extended. One future project 
is to look at HCN in more nearby galaxies, especially in individual regions forming super- 
star clusters, which may be the building blocks of star formation in starburst regions like 
ULIRGs and LIRGs (Ho 1997). ALMA will allow detailed study of HCN in other galaxies, 
including the J = 3 — > 2 transition. Another is to explore how far toward higher z these 
relations can be pushed, and to understand how the relations depend on metallicity and 
chemistry. Theoretical work by Lintott et al. (2005) suggest that HCN may trace very early 
star formation, depending on the nucleosynthetic yield of the earliest stars. 
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of Texas. YG acknowledges support of the NSF of China grants 10425313 & 10333060. 
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Fig. 1. 
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Fig. la 
points 



logL//j - logL^^^^_g correlation for Galactic and extragalactic sources, 
(above) indicates the linear least squares fit for Galactic cores (with Lm > 10^'^ L0 
above the dashed line) and for galaxies, separately. Fig. lb (down) shows a overall fit for 
both parts. The three isolated sohd squares are high z HCN 1-0 points from Solomon et al. 
(2003), Vanden Bout et al. (2004) and CariUi et al. (2005); they are not included in the fit 
because the sources are QSOs and the contribution from the AGN to Ljji is not yet clear. 
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Fig. 2. — Above: Correlation between Lm and Lir/Lhcni-o for galaxies (solid squares) and 
Galactic star forming cores (filled circles). Lir/Lhcni-o is constant for a large range of Ljr 
utill a turn off around Ljij=10^"^ Lq. Below: Lir per unit of molecular gas vs. molecular 
line luminosity. The star formation rate per amount of CO gas changes a lot from Galactic 
GMCs (the heavy line, with a dashed line boundary to indicate the variation), to galaxies 
(filled triangles) and high-z CO galaxies (hollow triangles). 



